TomograPy: A Fast, Instrument-Independent, Solar Tomography Software

Solar tomography has progressed rapidly in recent years thanks to the development of robust algorithms and the availability of more powerful computers. It can today provide crucial insights in solving issues related to the line-of-sight integration present in the data of solar imagers and coronagraphs. However, there remain challenges such as the increase of the available volume of data, the handling of the temporal evolution of the observed structures, and the heterogeneity of the data in multi-spacecraft studies. We present a generic software package that can perform fast tomographic inversions that scales linearly with the number of measurements, linearly with the length of the reconstruction cube (and not the number of voxels) and linearly with the number of cores and can use data from different sources and with a variety of physical models: TomograPy (http://nbarbey.github.com/TomograPy/), an open-source software freely available on the Python Package Index. For performance, TomograPy uses a parallelized-projection algorithm. It relies on the World Coordinate System standard to manage various data sources. A variety of inversion algorithms are provided to perform the tomographic-map estimation. A test suite is provided along with the code to ensure software quality. Since it makes use of the Siddon algorithm it is restricted to rectangular parallelepiped voxels but the spherical geometry of the corona can be handled through proper use of priors. We describe the main features of the code and show three practical examples of multi-spacecraft tomographic inversions using STEREO/EUVI and STEREO/COR1 data. Static and smoothly varying temporal evolution models are presented.Comment: 21 pages, 6 figures, 5 table

Automated Detection of Coronal Loops using a Wavelet Transform Modulus Maxima Method

We propose and test a wavelet transform modulus maxima method for the au- tomated detection and extraction of coronal loops in extreme ultraviolet images of the solar corona. This method decomposes an image into a number of size scales and tracks enhanced power along each ridge corresponding to a coronal loop at each scale. We compare the results across scales and suggest the optimum set of parameters to maximise completeness while minimising detection of noise. For a test coronal image, we compare the global statistics (e.g., number of loops at each length) to previous automated coronal-loop detection algorithms

Streamer Wave Events Observed in Solar Cycle 23

In this paper we conduct a data survey searching for well-defined streamer wave events observed by the Large Angle and Spectrometric Coronagraph (LASCO) on-board the Solar and Heliospheric Observatory (SOHO) throughout Solar Cycle 23. As a result, 8 candidate events are found and presented here. We compare different events and find that in most of them the driving CMEs ejecta are characterized by a high speed and a wide angular span, and the CME-streamer interactions occur generally along the flank of the streamer structure at an altitude no higher than the bottom of the field of view of LASCO C2. In addition, all front-side CMEs have accompanying flares. These common observational features shed light on the excitation conditions of streamer wave events. We also conduct a further analysis on one specific streamer wave event on 5 June 2003. The heliocentric distances of 4 wave troughs/crests at various exposure times are determined; they are then used to deduce the wave properties like period, wavelength, and phase speeds. It is found that both the period and wavelength increase gradually with the wave propagation along the streamer plasma sheet, and the phase speed of the preceding wave is generally faster than that of the trailing ones. The associated coronal seismological study yields the radial profiles of the Alfv\'en speed and magnetic field strength in the region surrounding the streamer plasma sheet. Both quantities show a general declining trend with time. This is interpreted as an observational manifestation of the recovering process of the CME-disturbed corona. It is also found that the Alfv\'enic critical point is at about 10 R$_\odot$ where the flow speed, which equals the Alfv\'en speed, is $\sim$ 200 km s$^{-1}$

The large longitudinal spread of solar energetic particles during the January 17, 2010 solar event

We investigate multi-spacecraft observations of the January 17, 2010 solar energetic particle event. Energetic electrons and protons have been observed over a remarkable large longitudinal range at the two STEREO spacecraft and SOHO suggesting a longitudinal spread of nearly 360 degrees at 1AU. The flaring active region, which was on the backside of the Sun as seen from Earth, was separated by more than 100 degrees in longitude from the magnetic footpoints of each of the three spacecraft. The event is characterized by strongly delayed energetic particle onsets with respect to the flare and only small or no anisotropies in the intensity measurements at all three locations. The presence of a coronal shock is evidenced by the observation of a type II radio burst from the Earth and STEREO B. In order to describe the observations in terms of particle transport in the interplanetary medium, including perpendicular diffusion, a 1D model describing the propagation along a magnetic field line (model 1) (Dr\"oge, 2003) and the 3D propagation model (model 2) by (Dr\"oge et al., 2010) including perpendicular diffusion in the interplanetary medium have been applied, respectively. While both models are capable of reproducing the observations, model 1 requires injection functions at the Sun of several hours. Model 2, which includes lateral transport in the solar wind, reveals high values for the ratio of perpendicular to parallel diffusion. Because we do not find evidence for unusual long injection functions at the Sun we favor a scenario with strong perpendicular transport in the interplanetary medium as explanation for the observations.Comment: The final publication is available at http://www.springerlink.co

Observations of quasi-periodic solar X-ray emission as a result of MHD oscillations in a system of multiple flare loops

We investigate the solar flare of 20 October 2002. The flare was accompanied by quasi-periodic pulsations (QPP) of both thermal and nonthermal hard X-ray emissions (HXR) observed by RHESSI in the 3-50 keV energy range. Analysis of the HXR time profiles in different energy channels made with the Lomb periodogram indicates two statistically significant time periods of about 16 and 36 seconds. The 36-second QPP were observed only in the nonthermal HXR emission in the impulsive phase of the flare. The 16-second QPP were more pronounced in the thermal HXR emission and were observed both in the impulsive and in the decay phases of the flare. Imaging analysis of the flare region, the determined time periods of the QPP and the estimated physical parameters of magnetic loops in the flare region allow us to interpret the observations as follows. 1) In the impulsive phase energy was released and electrons were accelerated by successive acts with the average time period of about 36 seconds in different parts of two spatially separated, but interacting loop systems of the flare region. 2) The 36-second periodicity of energy release could be caused by the action of fast MHD oscillations in the loops connecting these flaring sites. 3) During the first explosive acts of energy release the MHD oscillations (most probably the sausage mode) with time period of 16 seconds were excited in one system of the flare loops. 4) These oscillations were maintained by the subsequent explosive acts of energy release in the impulsive phase and were completely damped in the decay phase of the flare.Comment: 14 pages, 4 figure

Automated Coronal Hole Detection using Local Intensity Thresholding Techniques

We identify coronal holes using a histogram-based intensity thresholding technique and compare their properties to fast solar wind streams at three different points in the heliosphere. The thresholding technique was tested on EUV and X-ray images obtained using instruments onboard STEREO, SOHO and Hinode. The full-disk images were transformed into Lambert equal-area projection maps and partitioned into a series of overlapping sub-images from which local histograms were extracted. The histograms were used to determine the threshold for the low intensity regions, which were then classified as coronal holes or filaments using magnetograms from the SOHO/MDI. For all three instruments, the local thresholding algorithm was found to successfully determine coronal hole boundaries in a consistent manner. Coronal hole properties extracted using the segmentation algorithm were then compared with in situ measurements of the solar wind at 1 AU from ACE and STEREO. Our results indicate that flux tubes rooted in coronal holes expand super-radially within 1 AU and that larger (smaller) coronal holes result in longer (shorter) duration high-speed solar wind streams

An Extreme Solar Event of 20 January 2005: Properties of the Flare and the Origin of Energetic Particles

The extreme solar and SEP event of 20 January 2005 is analyzed from two perspectives. Firstly, we study features of the main phase of the flare, when the strongest emissions from microwaves up to 200 MeV gamma-rays were observed. Secondly, we relate our results to a long-standing controversy on the origin of SEPs arriving at Earth, i.e., acceleration in flares, or shocks ahead of CMEs. All emissions from microwaves up to 2.22 MeV line gamma-rays during the main flare phase originated within a compact structure located just above sunspot umbrae. A huge radio burst with a frequency maximum at 30 GHz was observed, indicating the presence of a large number of energetic electrons in strong magnetic fields. Thus, protons and electrons responsible for flare emissions during its main phase were accelerated within the magnetic field of the active region. The leading, impulsive parts of the GLE, and highest-energy gamma-rays identified with pi^0-decay emission, are similar and correspond in time. The origin of the pi^0-decay gamma-rays is argued to be the same as that of lower energy emissions. We estimate the sky-plane speed of the CME to be 2000-2600 km/s, i.e., high, but of the same order as preceding non-GLE-related CMEs from the same active region. Hence, the flare itself rather than the CME appears to determine the extreme nature of this event. We conclude that the acceleration, at least, to sub-relativistic energies, of electrons and protons, responsible for both the flare emissions and the leading spike of SEP/GLE by 07 UT, are likely to have occurred simultaneously within the flare region. We do not rule out a probable contribution from particles accelerated in the CME-driven shock for the leading GLE spike, which seemed to dominate later on.Comment: 34 pages, 14 Postscript figures. Solar Physics, accepted. A typo corrected. The original publication is available at http://www.springerlink.co

Physics of Solar Prominences: I - Spectral Diagnostics and Non-LTE Modelling

This review paper outlines background information and covers recent advances made via the analysis of spectra and images of prominence plasma and the increased sophistication of non-LTE (ie when there is a departure from Local Thermodynamic Equilibrium) radiative transfer models. We first describe the spectral inversion techniques that have been used to infer the plasma parameters important for the general properties of the prominence plasma in both its cool core and the hotter prominence-corona transition region. We also review studies devoted to the observation of bulk motions of the prominence plasma and to the determination of prominence mass. However, a simple inversion of spectroscopic data usually fails when the lines become optically thick at certain wavelengths. Therefore, complex non-LTE models become necessary. We thus present the basics of non-LTE radiative transfer theory and the associated multi-level radiative transfer problems. The main results of one- and two-dimensional models of the prominences and their fine-structures are presented. We then discuss the energy balance in various prominence models. Finally, we outline the outstanding observational and theoretical questions, and the directions for future progress in our understanding of solar prominences.Comment: 96 pages, 37 figures, Space Science Reviews. Some figures may have a better resolution in the published version. New version reflects minor changes brought after proof editin

Sources of SEP Acceleration during a Flare-CME Event

A high-speed halo-type coronal mass ejection (CME), associated with a GOES M4.6 soft X-ray flare in NOAA AR 0180 at S12W29 and an EIT wave and dimming, occurred on 9 November 2002. A complex radio event was observed during the same period. It included narrow-band fluctuations and frequency-drifting features in the metric wavelength range, type III burst groups at metric--hectometric wavelengths, and an interplanetary type II radio burst, which was visible in the dynamic radio spectrum below 14 MHz. To study the association of the recorded solar energetic particle (SEP) populations with the propagating CME and flaring, we perform a multi-wavelength analysis using radio spectral and imaging observations combined with white-light, EUV, hard X-ray, and magnetogram data. Velocity dispersion analysis of the particle distributions (SOHO and Wind in situ observations) provides estimates for the release times of electrons and protons. Our analysis indicates that proton acceleration was delayed compared to the electrons. The dynamics of the interplanetary type II burst identify the burst source as a bow shock created by the fast CME. The type III burst groups, with start times close to the estimated electron release times, trace electron beams travelling along open field lines into the interplanetary space. The type III bursts seem to encounter a steep density gradient as they overtake the type II shock front, resulting in an abrupt change in the frequency drift rate of the type III burst emission. Our study presents evidence in support of a scenario in which electrons are accelerated low in the corona behind the CME shock front, while protons are accelerated later, possibly at the CME bow shock high in the corona.Comment: Solar Physics, November 2007, in pres